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Resumen 

En este trabajo se plantea una nueva aplicación del modelo de las “ zetas ” con el fin de extender 
su utilidad a las correlaciones simultáneas de valores de la función de energía de Gibbs con la presión 
procedentes de equilibrios líquido-vapor isotérmicos, y los volúmenes de mezcla. E1 modelo de las 
“zetas ” se presenta como una relación indirecta con la concentración de la mezcla y con un parámetro k X) a 
través de la llamada fracción activa z. En el trabajo se analiza la efícacia del modelo en base a su 
extensión y a diferentes consideraciones realizadas sobre el parámetro k X] con el objetivo de conseguir la 
correlación más adecuada de los datos experimentales. Por último, se realiza una aplicación concreta 
sobre un conjunto de sistemas reales obtenidos en la bibliografía analizando los resultados de las 
diferentes propuestas e incluyendo un método para detenninar el parámetro /r¡j=\|/(/?). La función 
establecida para la función de energía de exceso de Gibbs tiene la forma genérica siguiente: 

^= z ( { ~ z )Yj b \ z ' siendo: B i=¿ 6 ¡j p 1 

KI i=0 

Palabras clave: Termodinámica; Correlación; Volumen exceso; Equilibrio líquido-vapor; Función de 
energía de Gibbs. 


Summary 

In this paper the authors plan a nevv application of the “ zetas ” model in order to extend its utility 
to the simultaneous correlations of the Gibbs energy values as a function of pressure, obtained from 
isotermal vapor-liquid equilibria, and the mixing volumes. The “ zetas ” model is presented as an indirect 
relationship with the mixing concentration and with a parameter k\ } through of the so called active fraction 
z. In this work the eficacy of the model considering its extension and difíerent forms to calculate the k i} 
parameter is analysed in order to get the more adequate correlation of experimental data. Lastly, a 
concrete application on a set of real systems extracted from bibliography is carried out analysing the 
results depending of different proposals and including a method to calculate the parameter £¡j=M/(/fí. The 
expression establisehd for excess Gibbs function presents the following generic form: 

^f= z ( l - z T£, B ¡ z ' where: B ¡=5Lhp ] 

Kl i=0 j=0 
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1. Introducción 

La energía de Gibbs es una función potencial termodinámica que se obtiene a partir de una doble 
transformación de la energía intema, en el siguiente orden. 

G = U-TS + PV = H-TS = A + PV (1) 

Debido a sus relaciones con otras cantidades termodinámicas, la función de energía de Gibbs se 
constituye como una de las fünciones potenciales termodinámicas más importantes. En ella se centran los 
estudios que venimos realizando en nuestro gmpo de investigación para su aplicación a casos de 
equilibrio de fases con soluciones, ya que, como veremos, las variables independientes de este potencial 
son medíbles fácilmente. Una de las ecuaciones derivadas de (1) sería la que relaciona a la fünción de 
Gibbs con la entalpía, o sea: 

G( p,T,N ) = H( p,T,N) + T 

identifícándose claramente el término del corchete con la fünción de entropía. 

Otra ecuación importante que surge de (1) sería la que permite el cálculo de la energía de Gibbs 
utilizando la relación que representa la dependencia con la presión de dicha fünción y de la energía de 
Helmholtz: 

G( p,T,N) = A( p,T,N ) + p 

la cual, de acuerdo con (1), produce la identifícación del corchete con V. Todas estas relaciones son 
aplicables, tal cual, a las correspondientes cantidades de exceso que surgen en los estudios rigurosos de 
disoluciones multicomponentes. 

La única propiedad extensiva que se recoge en las anteriores expresiones es el número de moles 
de las especies químicas presentes en el sistema en estudio. Con este trabajo y otros anteriores [1-3] se 
intenta plantear un procedimiento para determinar los valores de la fünción de la energía de Gibbs que 
surgen en los equilibrios entre fases de soluciones y otras cantidades termodinámicas, como son las 
entalpías y los volúmenes de mezcla. E1 método de cálculo con las entalpías fue presentado en un trabajo 
anterior [1], siendo objeto de este la utilización de datos de volúmenes y de equilibrios entre fases a 
temperatura constante. 

La importancia de la fünción de Gibbs en los estudios termodinámicos es clara ya que a partir de 
ella pueden calcularse los valores de las demás funciones termodinámicas y además, lo cual es posible 
porque es una fünción que se relaciona directamente con las magnitudes intensivas medibles por 
excelencia en el campo de la Termodinámica, como son la presión y la temperatura. Una vez planteada su 
relación con la entalpía y la entropía de exceso, a través de la propuesta de una expresión para la fünción 
de Gibbs con la concentración, se pretende introducir también la misma en las relaciones anteriores (1,3) 
con el fín de obtener unas expresiones directas con el volumen y la fünción de energía de Helmholtz A. 
Esta última fünción juega un papel destacado porque proporciona una conexión de la termodinámica con 
la mecánica estadística a través de la fúnción de partición. 

Los planteamientos de los casos comentados antes pueden llevarse a cabo de manera muy 
paralela, el primero se llevó a cabo en la publicación anterior [1] mientras que para el segundo caso se 
iniciará estableciendo la dependencia de la fünción del Gibbs con las magnitudes intensivas y de la/s 
cantidad/es de sustancia/s que intervienen en el sistema. Si un sistema fluido está constituido por varias 
sustancias y cada una de ellas interviene con las cantidades n u n 2 , .... , al ser la fünción de Gibbs una 
fünción homogénea de primer orden con respecto a estos parámetros extensivos puede escribirse: 

G -G(p,T,n v n 2 ,...) (4) 

Considerando una de las propiedades de las fünciones homogéneas y siendo n el número total de 
moles, con G=gn , la ecuación anterior también se escribe como: 


dG(p,T,N) 


d P 


(3) 


3G(pJ\Nl 

dT 


( 2 ) 
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g = g( p,T.x v x 2 ,...) 


(5) 


A partir de (4), o de (5), la variación que expenmenta la función de Gibbs con cada una de las 
variables, pucde cxpresarse matemáticamente por: 


dG = 


( dG} „ 

ícG] 

dp+y 


— dT + 



U T) pn 

U 'P ) 

1 L—i 

T.n ' 

l 5n J 


dn. 


( 6 ) 


p.T.n x 


expresión que puede aplicarse de la misma forma a las funciones de exceso. Después de operar 
adecuadamente y teniendo en cuenta la relación (1) y otros conceptos básicos de la Termociinámica cie 
Soiuciones , la (6) da lugar a una de las ecuaciones más importantes que se emplean para el estudio de los 
sistemas en equilibrio. 



dT+^—dp -^.v.c/ln/, =0 


(7) 


esta última ecuación es una versión de la anterior (6) pero con una caracterización más concreta de cada 
uno de los tres sumandos. Así, el primer sumando se asimila a la variación de la función de Gibbs con la 
tempcratura a través de un término que recoge la variación de entalpía de la solución; este aspecto fue 
analizado por Ortega y Espiau [1], lo mismo que el término correspondiente a la variación de G con la 
concentración, G=G(x), [2-4]. 

Aquí estableceremos un análisis termodinámico previo para plantear posterionnente una 
propuesta que recoja la \ ariación de la función de Gibbs con la presión, identificada en la ecuación (7) 
por los cambios de volumen de exceso de la solución considerada. La literatura recoge una densa 
información experimental sobre equilibrios entre fases, sobre todo de líquido-vapor obtenidos 
isotérmicamente, y que serán empleados en la valoración de la aplicación que se comenta. La 
organización del trabajo requiere que, en una primera parte, se plantee la relación que liga a las distintas 
fúnciones termodinámicas de mezcla (o de exceso) con la concentración v luego, en una segunda etapa, 

dirigir el artículo hacia la verifícación de la relación particular (<3g E / Sp\ x . que se corresponde con el 
segundo de los sumandos de la ecuación (6) o la (7). 


2. Representación empírica de la función de Gibbs-concentración-presión 

La representación de la fúnción de Gibbs con la concentración de un componente i en solución, 
g E =g E (.r¡), ha sido objeto de estudio y de aplicación por nuestra parte en trabajos anteriores [2-4]. La 
expresión genérica que se planteó inicialmente por Pacheco y Ortega [2] tiene la forma siguiente. 


p=2 


l|»2 -"'p 

¡,Í2—¡ P *CR*(n,p) 


( 8 ) 


ecuación que es una forma generalizada para una cantidad de exceso A f. Recordemos que con este 
modelo se establece una fúnción de exceso como suma de las contribuciones de los efectos de todas las 
posibles interacciones p-arias de las n fracciones activas , donde p= 2, 3, 4, etc. X representa el máximo 
orden de las interacciones moleculares que se consideren, mientras que los ¡ son los coeficientes 

particuiares de influencia de la combinación p-aria en la cantidad de exceso que se estudia: incluso, en 
algún caso uno o \ arios de esos coefícientes pueden ser nulos. La relación para \afraccióti activa r, de un 
determinado componente i se expresa también de manera genérica como: 

k x 

Z¡=-P- (i = 2,..., n) (9) 

x i + Í k n x > 

\=2 


siendo k ü un coeficiente a determinary al que se le ha dedicado un tratamiento diferenciado dependiendo 
de la propiedad. Se lle\ a a cabo una distinción considerando, bien un valor de A\ para correlacionar 


145 













volúmenes v e (jc), un valor de k h para entalpías /z E (x), que se relaciona con el anterior, ver [3], o bien un 
valor de k g para valores de la función de Gibbs. 

Para el caso de la solución más simple, como es la de una mezcla binaria, la ecuación (8) se 
simplifica en: 


= z, (l - Z| )k + A z, + ¿2 Z f ] 

JCJ 

donde ahora, z, referido a la sustancia “1”, viene a ser: z x =- 

X , +k 2i x 2 


( 10 ) 


( 11 ) 


Como se estableció en un trabajo anterior [1] existe una función para la energía de Gibbs de 
exceso que relaciona a la correspondiente función adimensional con la composición y con la temperatura, 
es decir, una relación del tipo g E =g E (T,x x ), que adoptó la forma más extendida siguiente: 


g E (x v T) 



+A i2 T + A i3 T 2 


( 12 ) 


Con la ecuación (12) y utilizando relaciones sencillas de la Termodinámica, se consiguieron expresiones 
para otras fúnciones termodinámicas, como h l y s E . 

De manera análoga, para plantear la dependencia con la presión, sería conveniente basamos en 
una de las clásicas relaciones de la Termodinámica por la que se relaciona a la función de Gibbs con la 
presión y con la temperatura. Dicha relación en su fomía diferencial establece que, para una sustancia 
pura que se constituye como un sistema termodinámico cerrado se cumple: dG = —SdT + Vdp , la que, 
particularizada a temperatura constante y para un fluido determinado con comportamiento ideal puede ser 
reescrita de la siguiente forma: 

dG = RT(dp/p ) = RT(d ln p) (13) 

la integración de esta sencilla ecuación, utilizada para introducir el concepto de fúgacidad de un fluido, 
podría ser: 

G = RT\np + A(T) (14) 


donde el término A(T) depende exclusivamente de la temperatura. De esta forma, la ecuación última 
agruparía a los dos sumandos originales, el de la presión y el de la temperatura. E1 desarrollo en serie de 
la función logarítmica ln p , en un entorno del punto ‘V’ sería: 

\np = \na + -(p-a)-^—(p-ay + ^—(p~ay -... (15) 

a 2 a 2 3 a 3 

de esta forma, la ecuación (14) de la función de Gibbs podría expresarse como una relación polinómica de 
la presión: 

G = RT\g^+g 2 p +g^p 1 +..]+A(T) (16) 

siendo los g¡ unos coeficientes que surgen del desarrollo polinómico de (15), por lo que, la dependencia 
de la función de Gibbs con la presión y la concentración, g E =g E (p, x), siguiendo un planteamiento análogo 
a aquel donde se estableció la dependencia con la temperatura [4], podria tener, considerando la (10), 

-- o^ =z » 0~*. ÍBo+B* +KA ] (17) 


pero donde ahora, de acuerdo a (16), la relación de B\ con p viene dada por: 

dB 

B=B n +B. 2 p+B a p 2 y —T=B a +2B B p 

dp 


(18) 
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( 19 ) 


- fk, + B oiP + B 0 ,P 2 )+ k i + B, 2 P + B nP' + ( fl :i + B nP + b zip' W 

Kl 

ahora la expresión anterior se puede plantear de manera más condensada como. 


RT 




Í(^P H )=l 


J 1 


( 20 ) 


La expresión (20) proporciona una buena aproximación de la función de Gibbs, si bien, por las 
consideraciones realizadas en esta propuesta, podrá tener unas limitaciones en cuanto al intenalo de 
presiones, pudiendo establecerse las restricciones correspondientes únicamente cuando se disponga de 
resultados procedentes de la aplicación a datos experimentales reales. Con la ecuación (20) y a partir de 
las relaciones (1) y/o a la (3) pueden obtenerse las relaciones para otras funciones termodinámicas, como 
v E . 

Así, de la primera derivada de (19) respecto a la presión, y teniendo en cuenta la identidad 
planteada en (3), considerando únicamente un polinomio en “r” de segundo orden para simplificación, 
resulta 


v E = 


dg 


E \ 


dp 


) T.x 


=z(\-z)j^ 


( dB , A 


dp 


z + 


dz Y dk 
Akdp 


(1-2 z)^B { z l +z(\-z)(B x +2 B 2 z) 


/L 


( 21 ) 


a partir de la relaciones escritas en (1) es claro que A = G — pV , la cual permite conseguir la siguiente 
ecuación 


£l-ll_ p 

RT RT 


Ís e /rt) 

cp 


resultando para la función de energía de Helmholtz. 


JT.x 


R 


=z(\-z)j 


B-p 

f dB > j 


(dtT 


d P J 

J \dk) 

\dp)_ 


{\-2z)jB,z'+z(\-zlB ]+ 2zB 2 ) 


( 22 ) 


En los modelos (21) y (22), el parámetro k influye más bien en la forma de la curva que en el tamaño de 
la ordenada, por lo que, en muchos casos cabe simplificar dichas ecuaciones, considerando que dicho 
parámetro es independiente de la presión, lo cual en este caso es bastante aceptable por tratarse de 
sustancias en estado líquido. De esta forma puede eliminarse el segundo de los sumandos de la ecuación 
(21) o de la (22), quedando una expresión de tipo polinómico sencilla en r, como la que se emplea 
usualmente en la correlación directa de las cantidades de exceso, \ er como ejemplo la ecuación (10). 


3. Aplicación del modelo propuesto 

Esta sección está dedicada al análisis de los resultados de ciertas aplicaciones sobre un conjunto 
de datos reales utilizando el modelo (19) para la fiinción de Gibbs, como ecuación de trabajo 
fúndamental, el cual puede escribirse en la forma condensada siguiente, 

,p\ 1-4*, .p)] ; (3) 

mientras, las funciones de exceso deri\ adas, como son el \ olumen y la función de Helmholtz, \ endrán 
dadas, respectivamente, por las ecuaciones (21) y (22). En cualquier caso, los coeficientes B¡{p) y sus 
derivadas matemáticas están expresadas en (18). 

Una vez presentado el formalismo termodinámico-matemático que se pretende en este trabajo, se 
está en condiciones de aplicar el modelo propuesto sobre aquellos sistemas reales para los que existan 
suficientes datos experimentales de equilibrio isotérmico líquido-\ apor y de volúmenes. Parece lógico 
que, tomando en consideración el modelo presentado anteriormente, pueden deducirse \ arios casos que 
implementaremos con el fín de encontrar la solución adecuada, todo ello dependiendo de la aplicación 
que se lle\ e a cabo. 
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La aplicación del procedimiento detallado en este trabajo y del modelo propuesto se realizará 
sobre un conjunto de sistemas binarios, principalmente con mezclas (alcanol+alcano), de los que se posee 
amplia información experimental. En nuestro laboratorio, a pesar de disponer de los equipos necesarios, 
no se han realizado determinaciones de equilibrio líquido-vapor (ELV) en condiciones isotérmicas, siendo 
esta la razón de acudir a bases de datos de la literatura intemacional con el fín de conseguir los datos 
precisos para valorar el método que se presenta. Resaltar que, a pesar de todo, no existe en la literatura un 
número considerable de sistemas que contenga la información necesaria para realizar una verifícación 
exhaustiva del procedimiento que se utiliza. 

E1 procedimiento a seguir parte de una recogida de datos experimentales, con un conjunto de 
valores (r, p , x, y) para cada uno de los estados de equilibrio, y de volúmenes de mezcla v E , en unas 
determinadas condiciones. E1 cálculo de los coeficientes de actividad y de los valores de la íúnción de 

Gibbs, g E / RT = ln/j , permite comenzar el proceso cuyo detalle se indica en el esquema de la 
fígura 1. Para un conjunto de mezclas binarias se llevó a cabo un ajuste simultáneo de las parejas de datos 
[x¡, v f ] y [ x \>/i(p)] utilizando un procedimiento de mínimos cuadrados, optimizando una fúnción objetivo 
FO defínida de la siguiente forma: 


ro=£ 


vL(X,’P)~ V Lexp 

2 

2 n 

. vv 

y í.cí ( p,’ x \,)~ /,j y ij..xp 

v E 

L ‘- ex p J 

2^2^ 

¡=1 j=i 



(24) 
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que representa la suma de los cuadrados de las diferencias entre las cantidades calculadas por el modelo y 
las obtenidas experimentalmente. Se observa que la ecuación (24) no contiene a los valores de la función 
dc Gibbs g L / RT , ya que estos se calculan a partir de los coeficientes de actividad y por tanto, 
estadísticamente no proporcionan una información independiente. Por el número de parámetros que se 
recogen en las ecuaciones mostradas y para evitar el problema de multiplicidad de raíces, se utilizó un 
algoritmo genético (AG) [5] que produce convergencia única para las soluciones. 

Con este procedimiento común de trabajo se han ideado diferentes propuestas o casos, cuyo 
análisis se pretende realizar a continuación, con el fin de concluir en la consecución de un método más 
óptimo, tanto matemática como termodinámicamente. Los diferentes casos que pueden presentarse así 
como las correspondientes aplicaciones de cada uno de ellos, se detallan en las siguientes secciones. 


3.1. Empleo de un modelo simplificado 

Este caso se refiere a la aplicación del modelo propuesto teniendo en cuenta únicamente un 
número reducido de coeficientes que surge de truncar el número de términos de la primera de las 
expresiones (18) en solo dos. Pueden plantearse dos casos que denominaremos 1 y 2, dependiendo de si el 
parámetro £¡j es, respectivamente, dependiente o no de la presión, tal como se indica en el esquema de la 
figura 1. E1 denominado caso 1 se comentará en la sección siguiente. En el segundo de los casos, el valor 
de /:¡j, identificado por k v para los v 1 , o con k g para la g l /RT, se considerarán idénticos e independientes 
de la presión , y por tanto constantes en la correlación planteada. Es decir, los valores de dicho parámetro 
se calcularán de acuerdo con la ecuación (A.l) del apéndice, quedando clara la trayectoria seguida en el 
diagrama de flujo indicado en el esquema de la figura 1, modelo simplificado->/r v ^(p(p)-+e cuación 

(A.l), vía 2a. Otra posibilidad para este caso, que corresponde a la vía que aparece en la otra rama del 
esquema, la 2b, sería la de conseguir los valores del parámetro k s .~k g en el mismo procedimiento de 
regresión, como un parámetro más, al correlacionar conjuntamente los datos de ELV isotérmicos y los v E . 
Lógicamente, en estos casos se cumplirá que ( ák- / áp) = 0 , por lo que las ecuaciones (18), (19), (21) y 
(22) pueden adquirir, respectivamente, las formas simplificadas siguiente: 


B \ ~ B ¡o +B ¡\ P 


y 


- L = B, , entonces, 

d P " 


(25) 


E 2 

—• = + B m p)+{B w + , p)z +(£, 0 + B 2] p)z ]= z(\ — z)^(5 ¡0 + B it p)z (26) 


RT 

v E 

Jt' 


dj^RT) 

dp 


=z(l-zfB 0i +B n z+B 2i z 2 ]= z(\-z)j^B n z' ¡ 


(27) 


-Ir.jr 




B-p 


z' =z(\-zÍB w +B m z+B 10 z 2 ) 


(28) 


Los resultados de las dos aplicaciones a las que se ha hecho referencia sobre el conjunto de datos 
elegidos, se recogen en las tablas l(a) y l(b). Así, la primera de las tablas muestra los coeficientes 
obtenidos para un conjunto de sistemas, donde el parámetro k x =kz se consideró fíjo en el proceso de 
correlación, calculado por la ecuación (A.l), y aparece en la primera columna numérica de la tabla l(a), 
acorde a la trayectoria 2a del esquema de la fígura 1. Finalmente, la aplicación del procedimiento de 
regresión, tal como se indicó en la sección anterior, permitió completar las demás columnas de la 
mencionada tabla, incluyéndose los parámetros de bondad de ajuste que se indican. 

En la tabla l(b) se recogen los valores obtenidos para el segundo caso, la trayectoria 2b. Como se 
indicó, ahora el parámetro k s ^k g se somete al procedimiento de regresión del AG siendo uno más a 
determinar en el proceso de correlación indicado en la § 3, obteniéndose fmalmente un conjunto único de 
parámetros para las ecuaciones (26) y (27) que permiten su adecuada representación. 
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Tabla 1 (a). Valores de Á¡¡ y B \¡ obtenidos para el modelo (26), en la correlación conjunta de los ELV 
isotérinico y los v' a la temperatura indicada (7), de diferentes mezclas binarias, siendo s( y) la desviación 
estándar de la cantidad “y”. Aquí Q=(g L /RT) y k g = k K , indicado como k l} , se obtiene mediante (A.l) y es 
fijo en el procedimiento de ajuste. 


Mezcla 

k 

Bqo 

Z?oi x 10 3 B\q Z?]i x 10 3 

B 20 

£ 2 i*10 3 

s(Q) 

s(?i) 

X(V L ) 

^ietanol + 

v 2 hexano (298.15 ) 6 ' 

2.242 

5.360 

1.746 -9.594 -3.737 

5.852 

2.607 

0.033 

1.138 

0.036 

jr 2 hexano (318.15) 74 

2.256 

6.036 

2.410 -11.44 -6.004 

6.626 

4.474 

0.075 

3.270 

0.059 

x 2 heptano (298.15 ) 6 ‘ 8 

2.514 

4.584 

2.056 -6.899 -4.141 

3.218 

2.702 

0.058 

1.057 

0.048 

x 2 heptano (313.15) 9 ' 10 

2.519 

5.076 

2.244 -9.366 -4.760 

5.945 

3.343 

0.029 

0.768 

0.042 

vipropan-l-ol + 

jiheptano (298.15)"'" 

1.963 

4.242 

1.279 -6.174 -2.414 

3.035 

1.366 

0.026 

0.473 

0.025 

x 2 nonano (298.15) 12 ' 3 

2.391 

4.414 

1.687 -7.143 -3.398 

3.978 

2.127 

0.021 

0.315 

0.040 

xibutan-l-ol + 

x 2 hexano (298.15) 1415 

1.431 

2.982 

0.585 -3.108 -1.688 

1.112 

1.008 

0.023 

0.806 

0.021 

v 2 heptano (313.15) 9 ' 16 

1.635 

4.126 

1.360 -6.615 -3.337 

3.808 

2.291 

0.030 

1.214 

0.027 

jcjisobutil alcohol + 

x 2 hexano (298.15) 1417 

1.416 

0.926 

1.027 0.914 -2.204 

-1.066 

1.247 

0.069 

0.497 

0.012 

jc 2 heptano (298.15) 12 ' 18 

1.588 

3.448 

1.435 -4.161 -2.734 

1.859 

1.522 

0.016 

0.213 

0.024 

x 2 heptano (313.15) 16 

1.594 

3.323 

1.692 -3.871 -3.061 

1.627 

1.637 

0.022 

0.971 

0.025 

A'idimetil carbonato + 

jc 2 hexano (298.15) 19 ' 20 

1.554 

2.959 

3.253 -3.194 -6.128 

1.784 

3.660 

0.005 

0.067 

0.021 

j 2 octano (298.15) 19 20 

1.931 

3.263 

4.606 -4.024 -7.936 

2.249 

4.482 

0.009 

0.082 

0.024 

x 2 decano (298.15) 19 ' 2 

2.313 

3.934 

5.310 -5.975 -9.823 

3.853 

6.405 

0.031 

0.200 

0.045 


Tabla 1 (b). Valores de £¡¡ y B t 

, obtenidos para el modelo (26), 

en la correlación conjunta de los ELV 

isotérmico y los v a 

la temperatura indicada (7), de diferentes mezclas binarias, 

siendo 

s(y) la 

1 desviaeión estándai de la eantidad ~y’ 

‘. Aquí y-(g F //?f) > A g 

- A v , indicado coino A, 

se obtiene por 

optimización de (24) en 

i el procedimiento de ajuste. 






Mezcla 

k 

Bqq 

Z?oi x 10 3 R\o ^ii x 10 3 

#20 

S 21 xl0 3 

s(Q) 

■s(j'i) 

í(v e ) 

Jietanol + 

1 x 2 hexano (298.15) 6 ' 4 

0.743 

2.416 

0.969 -1.549 -1.266 

2.308 

1.289 

0.020 

0.190 

0.008 

j 2 hexano (318.15) 7 ' 4 

0.641 

2.159 

1.180 -1.845 -1.303 

3.151 

1.136 

0.004 

0.338 

0.016 

jc 2 heptano (298.15) 6 ' 8 

0.785 

1.946 

1.112 -1.555 -2.462 

3.831 

3.233 

0.011 

0.391 

0.024 

jc 2 heptano (313.15) 9 ' 10 

0.863 

2.453 

1.399 -1.285 -1.976 

1.623 

1.849 

0.001 

0.034 

0.004 

A'ipropan-l-ol + 
jc 2 heptano (298.15) 118 

0.651 

2.108 

0.924 -1.216 -1.387 

2.237 

1.091 

0.003 

0.101 

0.006 

jc 2 nonano (298.15) 12 * 13 

0.922 

2.421 

1.145 -1.457 -1.750 

1.739 

1.444 

0.003 

0.034 

0.011 

jcjbutan-l-ol + 
jc 2 hexano (298.15) 14,15 

0.556 

1.772 

0.447 -0.762 -0.309 

1.421 

-0.447 

0.012 

0.227 

0.013 

jc 2 heptano (313.15) 916 
jci isobutil alcohol + 

0.617 

1.939 

0.771 -0.856 -0.774 

1.657 

0.271 

0.016 

0.109 

0.013 

jc 2 hexano (298.15) 1417 

0.441 

0.511 

-0.096 1.156 2.023 

0.923 

-1.892 

0.015 

0.114 

0.032 

j 2 heptano (298.15) 12 ' 18 

0.605 

1.835 

0.661 -0.757 0.027 

1.815 

-0.241 

0.003 

0.076 

0.015 

jc 2 heptano (313.15) 16 
jcidimetil carbonato + 

0.618 

1.611 

0.800 -0.337 0.063 

1.449 

-0.291 

0.007 

0.568 

0.016 

jc 2 hexano (298.15) 19 ' 20 

0.889 

1.809 

2.124 -0.336 -1.812 

0.870 

0.610 

0.001 

0.014 

0.005 

x 2 octano (298.15) 19 ’ 20 

0.965 

1.867 

2.533 -0.417 -0.690 

0.962 

-0.148 

0.001 

0.004 

0.005 

jc 2 decano (298.15) 19 ' 20 

1.063 

2.374 

3.093 -1.611 -2.579 

1.757 

2.219 

0.001 

0.006 

0.010 


La comparación de los dos casos mencionados en esta sección, permiten destacar importantes 
diferencias que pueden resumirse. E1 segundo de los casos, vía 2b, con el valor del parámetro k\¡ incluido 
en el mismo procedimiento de correlación y optimizado por el, da lugar a correlaciones sensiblemente 
mejores de todas las cantidades implicadas, como son los coeficientes de actividad y;, la íunción 
adimensional de Gibbs Q=g E /RT , y los volúmenes de exceso v E , comprobándose este hecho para todos los 
casos estudiados y especialmente para la mezcla (etanol+hexano) a 318,15 K, remarcado en las tablas. 
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3.2. Empleo de un modelo simplificado con dependiente de p 

Se trata en este caso utilizar el mismo modelo anterior (26), pero distinguiendo entre el valor de 
utilizado para los volúmenes y para los datos de ELV, es decir, k s ^ k & . Ahora el \ alor de k s se hace 

depender de la presión de acuerdo con la ecuación (A.4) del apéndice, todo ello si se consiguen \ alores de 
volúmenes a diferentes presiones (bastante poco usual), mientras que el valor de k¿, indicado para el 
ajuste de la función de Gibbs, se conseguirá de acuerdo al proceso indicado en el apcndice, estableciendo 
previamente unos \ alores iniciales conseguidos mediante (A.l) el de A\, y el de utilizando datos de 
ELV isobárico a /? 0 = Gte. Luego, las \ ariaciones de este parámetro con p se plantean, respectivamente, 
según (A.4) y (A.7). La aplicación del algoritmo a los sistemas elegidos produce los resultados de la tabla 
l(c). Si se comparan los resultados de esta tabla con los de las anteriores, se aprecia que para este caso las 
correlaciones no demuestran una mejoría significativa en la utilización del modelo simplificado con los 
parámetros k l} dependiente de p (lo cual parece bastante lógico), notándose una mejoría relati\ a solo en 
algunos casos para la fúnción adimensional de Gibbs y para los coeficientes de acti\ idad, lo cual también 
era de esperar por ser ahora el parámetro k § dependiente de cada estado de equilibrio. Por ello, el resumen 
de esta aplicación sería clasificarlo como un caso intermedio entre los dos anteriores. 


Tabla 1 (c). Valores de k x } y #¡j obtenidos para el modelo (26), en la correlación conjunta de los ELV 
isotérmico y los v E a la temperatura indicada (T), de diferentes mezclas binarias, siendo 5(y) la desviación 
estándar de la cantidad “y”. Aquí Q=(g E /R7 ). Los \ alores de k y y k % son los que se plantean como valores 
iniciales para su empleo en (A.4) y (A.7). 


Mezcla 

ko(p 0 ) 

A\(po 

) 

£oo 

#01*10 #i 0 #,,xlO' 

#20 

B : ,X10 3 

s(Q) 

■S(}’¡) 

s(l’ fc ) 

jc^etanol + 










jc 2 hexano 

1.767 

2.242 

4.252 

1.765 -6.032 -3.828 

3.557 

2.698 

0.02 

0.732 

0.036 

(298.15) 64 







3 



jc 2 hexano (318.15) 7 ' 4 

1.767 

2.256 

4.701 

2.346 -7.458 -5.637 

4.304 

4.058 

0.04 

0 

2.229 

0.058 

x ; heptano (298.15) 6 " 

1.973 

2.514 

4.416 

2.030 -6.495 -4.484 

3.766 

3.334 

0.03 

5 

0.01 

0.567 

0.041 

jc 2 heptano 

1.973 

2.519 

4.347 

2.248 -6.606 -4.781 

3.912 

3.365 

0.491 

0.042 

(313.15) 910 
jc^propan-l-ol + 







4 



jc 2 heptano 

1.615 

1.963 

3.763 

1.361 -5.072 -2.933 

2.935 

1.941 

0.01 

0.267 

0.021 

(298.15)" 8 







4 



x ; nonano 

1.967 

2.391 

4.065 

1.744 -6.123 -4.049 

3.770 

3.043 

0.01 

0.200 

0.033 

C298.15) 1213 







6 



jcjbutan-l-ol + 
.v 2 hexano 

1.219 

1.431 

2.954 

0.743 -3.982 -2.227 

2.707 

1.463 

0.00 

0.337 

0.016 

(298.15) 1415 







8 



jc 2 heptano 

1.370 

1.612 

3.584 

1.351 -4.899 -3.288 

2.823 

2.243 

0.02 

0.873 

0.026 

(j1 j>. 15) 







3 



jcjisobutil alcohol + 
jc 2 hexano 

1.212 

1.416 

1.764 

1.065 -0.828 -2.432 

0.153 

1.497 

0.00 

0.110 

0.010 

(298.15) 1417 







8 



jc 2 heptano 

1.358 

1.588 

3.218 

1.435 -3.763 -2.960 

2.152 

1.864 

0.00 

0.103 

0.020 

(298.15) 121 * 







8 



x ; heptano (313.15) 16 

1.361 

1.594 

3.332 

1.724 -4.210 -3.381 

2.344 

2.055 

0.01 

Q 

0.461 

0.021 

v^dimetil carbonato+ 







O 



.v 2 hexano 

1.264 

1.554 

2.466 

3.253 -1.727 -6.122 

1.014 

3.652 

0.00 

0.029 

0.021 

(298.15) 1920 







2 



xioctano 

1.720 

1.931 

2.992 

4.611 -3.140 -7.961 

1.703 

4.507 

0.00 

0.058 

0.024 

(298.15) 1920 







6 



v 2 decano 

2.163 

2.313 

3.808 

5.322 -5.588 -9.881 

3.637 

6.461 

0.02 

0.169 

0.045 

(298.15)' 020 







7 
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3.3. Empleo de un modelo extendido con k V] independiente de p 

Calificamos como modelo extendido a aquel que se plantea teniendo en cuenta un número mayor 
de términos en p para los coeficientes B x (p\ ver ecuaciones (18) a (22). De esta forma, dichos 
coeficientes, sus derivadas y el modelo para cada una de las cantidades de mezcla, adoptan las 
expresiones concretas que se escriben a continuación. Para esta primera aplicación también se han 
considerado independientes de la presión los parámetros k s . y k g , pero haciéndolos idénticos . y anulándose 
las correspondientes derivadas \dk Vj /dp). Aquí es aplicable todo lo que se detalló en la § 3.1, 

incluyendo lo referente a las dos vías de cálculo que se proponen en el diagrama de flujo de la fígura 1. 
Las ecuaciones para este caso son: 


B \ ~B\o + Bi\P + B-aP 2 +B¡iP 3 


dB, 


dT 


r-^i + 2#i 2 P + ^ B aP~ 


r E 

RT 

,E 


= z(l - z)£ (B i 0 +B a p+B n p 2 +B n p 2 )z' 


1=0 

2 


^-=z(l-z)$>, + 2B aP+ 3B n p 2 )z' 

A I ;=n 




RT 


■A l ~ z )i( B \o ~ B aP 2 +2 b , 3 p' )z' 


(29) 

(30) 


(31) 


(32) 


Los resultados obtenidos en la aplicación de esta versión más “extendida” del modelo sobre el mismo 
conjunto de sistemas que se han empleado en los casos anteriores se muestran en las tablas 2(a) y 2(b). 
Un planteamiento paralelo al de la § 3.1, tal como se recoge en el esquema de la fígura 1, sería el 
siguiente: un primer caso, vía 2, modelo extendido—» k g =k y ^ (p(p) ->ecuación (A.l), recogiéndose 

los resultados de esta aplicación en la tabla 2(a); puede verse que los valores de k x coinciden exactamente 
a los de la tabla l(a). E1 segundo caso, vía 2b, corresponde a la trayectoria del diagrama: modelo 
extendido ->k g =k y ^(p(p)->v alor constante, pero obtenido en el ajuste simultáneo con todos los 

datos empleando el AG; los resultados finales aparecen en la tabla 2(b). 

Los resultados obtenidos con la versión extendida, tablas 2(a) y 2(b), ofrecen una cierta mejoría 
en relación a los mostrados en las tablas correspondientes a la versión simplifícada, l(a) y l(b), 
especialmente para aquellos casos considerados como “más difíciles”, como son los del (butan-1- 
ol+hexano) y (2-metilpropan-l-ol+hexano o +heptano). Todo ello como consecuencia de aumentar el 
número de coeficientes del modelo. A1 igual que en el modelo simplifícado también aquí se pone de 
manifiesto la mejora que supone optimizar la correlación, calculando el parámetro k X] en el procedimiento 
de ajuste. 


3.4. Empleo de un modelo extendido con k V] dependiente de p 

En este caso se utiliza el modelo (30) para correlacionar los valores experimentales de los 
sistemas elegidos y el valor de k v con una dependencia de la presión según (A.4) y fijando un valor inicial 
de k x en condiciones de p=p 0 . Mientras, el valor de k % se hace variar con la presión según la (A.7) a partir 
de un valor inicial obtenido en la correlación del mismo sistema binario en condiciones isobáricas de 
p=Po- Los ajustes obtenidos con este procedimiento aparecen en la tabla 2(c). Según el diagrama de flujo 
de la fígura 1, este caso corresponde a la trayectoria: modelo extendido-> k V} =(p(p) —>ecuaciones (A.4) 
y (A.7). Los resultados de la correlación indican que el aumento del número de coefícientes no 
proporciona una mejoría relevante con respecto a los valores obtenidos en el caso homólogo de la tabla 
l(c), únicamente algunos casos puntuales presentan una desviación menor, tanto en los volúmenes de 
exceso como en la función de Gibbs, que no compensa una mayor difícultad para el modelo. 
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Tabla 2 (a). Valores de k y B l} obtenidos para el modelo (29), en la correlación conjunta de los ELV isotérmico v los v a la 
temperatura indicada (T), de diferentes mezclas binarias, siendo s(v) la desviación estándar de la cantidady. Aquí Q={¿ RT) 
v k e = k„ indicado como k , se obtiene mediante (A.1.1) y es fijo en el procedimiento de ajuste. 


#:¡ 

Bci j B | ( B' ¡ 

#Oi 


B ,, 

B S 


.i'ietanol + .v : hexano 

j,etanol + .v : hexano 

x,etanol+x : heptano 

x,etanol+x : heptano 


(298.15) 64 

(318.15) 4 

(298.15r 


(313.15) 




6.415 - 6.232 



5.131 -9.442 

5.958 


5.366 -9.671 5.867 


5.244 -9.184 

5.726 



/=o 


11.755 






5.124* 1.691* 1.690* 

. 3 . 44 * 4.444 1.689* 

4.547 -8.226 

6.254* 

1.250* 9.020* 

1.690* 

/= i 

10 J I0' 5 I0' 3 

I0 4 *10' 3 I0' 2 

*10' 3 *10' 3 

10* 3 

10* 4 10* 4 

10' 3 


-6.851 -6.262 2.983* 

-1.042 -6.225 1.537* 

7.001* -6.395 

-5.550 

-9.189 -1.849 

3.138* 

y-2 

xlO -6 xlC 6 I0' 5 

x|0 4 x|0' 5 I0 4 

10* 6 40* 6 

xlO -6 

40* 6 *10' 5 

10* 5 


8.582* -1.381 -1.635 

7.741* 3.654* -1.424 

1.676* 1.630* 

-5.773 

1.296 * 4.933 

-1.515 

7=3 

I0' 8 x|0' 7 x|0' 7 

I0' 7 I0' 7 xio- 6 

10' 7 10' 7 

x|0' 8 

10* 7 *10' 9 

*io* 7 

í(0 

0.034 

0.070 

0.053 


0.029 


i(/¡) 

1.132 

3.111 

0.809 


0.765 


s(v l ) 

0.036 

0.058 

0.041 


0.042 



.r | propan-1 -ol+x : heptano 

.v, propan-1 -ol+.vmonano 

x,butan-1 -ol+x-hexano 

x,butan-l -okv : heptano 


(298.15)" 8 

(298.l5)' lb 

(298.15) 1415 

(313.15) 9 ' 6 


7=0 

4.380 -7.125 4.567 

4.367 -7.304 4.535 

3.511 -5.174 

3.065 

4.095 -6.568 

3.786 


-4.533 -1.665 -1.681 

-6.160 6.408* 6.889* 

-7.480 -5.486 

7.823* 

1.150* 1.691* 

1.690* 

7=1 

xio -1 x 10 2 x|0' : 

x|0 4 I0' 3 I0' 3 

*10* 3 *10* 3 

10* 3 

10* 3 10* 3 

10* 3 


-3.346 -1.555 -1.551 

3.496* -2.762 1.075* 

1.496* 6.644* 

-8.155 

9.634* -6.650 

2.822* 

7=2 

x|0' 5 *10 4 xio 4 

I0' 5 x|0' 5 I0- 6 

10* 5 10* 8 

40* 6 

10* 8 xlO* 6 

10* 5 


2.792* 1.468* 1.630* 

-1.533 -1.581 -1.316 

1.690* 1.036* 

-1.511 

5.791* -1.173 

-1.684 

/=3 

10 7 iO 6 10* 6 

*10' 7 *10' 7 *10' 7 

10* 7 10* 7 

*10* 7 

10* 9 *10* 7 

*10 7 

s(P) 

0.024 

0.024 

0.013 


0.030 



0.404 1 

0.291 

0.596 


1.188 



0.021 

0.033 

0.016 


0.026 



.v,isobutil alcohol+v : he.\ano 

A,isobutil alcohol+v : heptano x,isobutil alcohoh.v ; 

; heptano 

r,dimet¡lcarbonatorc 

: he.\ano 


(298.I5) 1417 

(298.15) 8 

(313.1 5f 


(298.15) 19 - 0 ' 


7=0 

2.217 -2.325 1.056 

3.685 -5.144 3.048 

3.780 -6.102 

3.403 

2.967 -3.112 

1.758 


1.367* -7.490 -8.450 

-8.200 -7.282 -8.403 

5.443* 1.690* 

1.492* 

2.381* -8.076 

1.390* 

7=1 

10' 3 *10‘ 5 *10' 3 

*10' 3 *10' 3 * 10' 3 

10* 3 10* 2 

10‘ 2 

10* 3 *10* 3 

10' 3 


2.418* -8.988 2.556* 

3.627* -1.765 2.504* 

-1.970 -1.024 

-3.760 

1.631* -6.467 

9.773* 

7=2 

10' 5 *10' 7 10' 5 

10* 5 ^IO* 6 10' 5 

*10' 5 *10 4 

*10* 5 

10* 5 *10* 5 

10* 5 


-1.687 -7.166 1.560* 

7.430* 1.515* 0.169* 

9.246* 1.341* 

1.691* 

-7.896 4.885* 

-5.691 

7=3 

*10' 7 *10' 8 10' 7 

10' 8 10' 7 10' 7 

10* 9 10* 8 

10* 7 

*10* 8 10‘ 7 

*10* 7 

í(0 

0.007 

. 0.013 

0.027 


0.005 


í(J'í) 

0.092 

0.163 

0.624 


0.067 


s(v E ) 

0.010 

0.020 

0.022 


0.021 



.V|dimetilcarbonato+ 

* x,dimetilcarbonato+ 






x^octano 

x-decano 






(298.I5) 19 - 0 

(298.I5) 19 " 0 





7=0 

3.284 4.041 2.268 

3.947 -6.046 3.844 






7.353* -2.469 -8.963 

3.564* 3.854* 8.451* 





7=1 

I0" 4 *10' 3 40 4 

10' 3 10* 4 10' 3 






3.223* -3.083 4.993* 

7.439* -3.552 1.171* 





7=2 

10* 6 *10' 5 10' 5 

10* 6 *10' 5 10* 5 






1.048* 2.367* -1.524 

8.391* -1.007 -1.406 





7=3 

10' 7 10' 8 *10' 7 

10' 9 *10' 7 *10' 7 





m 

0.009 

0.031 





s(’/¡) 

0.081 

0.199 





s(v h ) 

0.024 

0.045 
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Tabla 2 (b). Valores de k y B i} obtenidos para el modelo (29), en la correlación conjunta de los ELV isotérmico y los v h a la 
temperatura indicada (T), de diferentes mezclas binarias, siendo s(y) la desviación estándar de la cantidad y. Aquí 
Q^i^lRT) y k g = k v , indicado como k, se obtiene mediante (A.1.4) y es fijo en el procedimiento de ajuste. Aquí Q=(¿'/RT) 
yk a = k^ , indicado como k, se obtiene por optimización de (23) en el procedimiento de ajuste. 


B()\ 

«n 

§2i 

B o\ 

8n 


Bo\ 

8n 


Bq\ 

§ii 

Bii 

x^etanol+xnhexano 

(298.15) 6 ' 4 

jC]etanol+x 2 hexano 

(318.15) 7 ' 4 

jcietanol+jc^heptano 

(298.15) 6 ® 

X] etanol+x 2 heptano 
(313.15) 9 ' 10 

7=0 

2.419 

-1.538 

2.345 

2.121 

-1.122 

3.289 

2.335 

-1.397 

2.496 

2.427 

-1.169 

1.660 


-1.348 

8.467x 

4.900x 

-2.763 

-5.703 

-6.760 

-1.395 

4.808x 

2.620x 

2.257x 

-1.415 

3.916x 

7=1 

XlO' 3 

10 4 

10' 4 

XlO' 3 

XlO' 3 

XlO' 3 

xl0‘ 3 

10* 4 

10 4 

10 4 

xl0‘ 3 

10' 4 


4.132x 

5.796x 

5.892x 

1.414x 

2.435x 

1.393x 

-6.718 

-7.826 

-1.271 

-1.267 

-1214 

-1.690 

7=2 

10' 6 

10' 6 

10' 6 

10' 5 

10' 7 

10* 5 

xlO' 6 

xlO' 6 

xlO* 5 

xlO* 5 

xlO" 5 

xlO" 5 


4.786x 

-1.068 

-1.200 

3.355x 

1.417x 

1.690x 

1.270x 

-2.476 

1.439x 

1.209x 

6.309x 

1.584x 

7=3 

10 ® 

XlO' 7 

Xl0' 8 

10' 8 

io - 7 

10' 7 

10‘ 7 

xlO' 8 

10’ 7 

10‘ 7 

10‘ 8 

10' 7 

m 

0.020 


k = 
0.729 

0.004 


k = 
0.560 

0.015 


k = 
0.705 

0.001 


k = 
0.840 

s(y¡) 

0.184 



0.234 



0.102 



0.031 




0.008 



0.014 



0.010 



0.004 




jc, propan-1 -ol+jc 2 heptano 
(298.15) 11 - 8 

jcipropan-1 -ol+jc 2 nonano 
(298.15) 12 ' 13 

jcjbutan-1 -ol+jc 2 hexano 
(298.15) 14 ' 15 

xibutan-1 -ol+x 2 heptano 
(313.15) 9 ' 16 

j =o 

2.039 

-0.785 

1.717 

2.556 

-1.705 

1.725 

2.089 

-1.816 

2.186 

1.904 

-0.784 

1.698 


-5.507 

4.889x 

-1.690 

1 .202x 

-1.689 

-1.252 

-2.679 

-6.717 

1.213x 

-4.925 

-9.826 

-1.646 

7=1 

XlO' 5 

10' 4 

XlO' 3 

10' 4 

XlO' 3 

XlO' 3 

xlO' 6 

xlO 4 

10' 3 

xl0‘ 4 

xlO" 4 

xlO' 3 


-6.516 

-6.801 

-9.255 

-1.889 

-4.436 

1.242 

-1.676 

-1.689 

-1.658 

3.173x 

2.380x 

2.982x 

7=2 

Xl0' 7 

Xl0' 6 

xlO' 6 

XlO 4 

X|0 6 

XlO' 6 

xlO' 5 

xlO' 5 

xlO' 5 

10' 5 

10' 5 

10' 5 


3.355x 

-7.215 

1.438x 

4.485x 

3.056x 

7.614x 

1.260x 

1.128x 

6.573x 

-1.683 

-1.179 

-1.352 

7=3 

10' 8 

XlO' 9 

10' 7 

10' 8 

10' 8 

10' 8 

10' 7 

10' 7 

10' 8 

xlO' 7 

xlO' 7 

xlO' 7 

í(5) 



k = 



k = 



k = 

0.016 


k = 

0.004 


0.672 

0.003 


0.980 

0.003 


0.653 



0.601 

síyi) 

0.025 



0.018 



0.040 



0.105 



s(v h ) 

0.006 



0.010 



0.011 



0.013 




X|isobutil alcohol+x 2 hexano 
(298.15) 141S 

Xjisobutil alcohol+x 2 heptano X]isobutil alcohol+x 2 heptano 
(298.15) 121 ® (313.15) 16 " 

x 1 dimetilcarbonato+ 
x 2 hexano (298.15) 19 ' 20 

7=0 

1.500 

-0.147 

-0.108 

2.162 

-1.191 

1.591 

1.948 

-0.928 

1.545 

1.782 

-0.300 

0.921 


1.633x 

-1.691 

-1.690 

-1.321 

-1.352 

1.160x 

1.690x 

-2.707 

-4.668 

1 .200x 

3.565x 

-1.589 

7=1 

10’ 3 

xlO’ 3 

xlO' 3 

xlO' 3 

xlO' 3 

10* 3 

10* 3 

xlO' 4 

xlO* 4 

10‘ 3 

10* 4 

xlO' 3 


1.424x 

1.571 x 

-1.225 

-7.132 

-1.405 

-1.637 

1.426x 

-5.307 

-7.400 

1.322x 

-8.062 

-1.492 

7=2 

10‘ 5 

10‘ 5 

xlO' 5 

xio- 6 

xlO* 5 

xlO" 5 

10‘ 5 

xl0‘ 6 

xl0‘ 7 

10‘ 5 

xlO' 6 

xlO' 5 


-1.197 

-9.235 

1.542x 

1.175 x 

1.125x 

8.104x 

-1.170 

2.246x 

2.955x 

-5.824 

-1.372 

1.679x 

7=3 

xlO* 7 

xlO' 8 

ÍO' 7 

io - 7 

10' 7 

10' 8 

xlO' 7 

ÍO’ 8 

10' 8 

xlO' 8 

xlO' 8 

io - 7 

s(Q) 

0.008 


k = 



k = 



k = 



k = 



0.967 

0.002 


0.712 

0.012 


0.666 

0.0010 


0.869 

í(}’¡) 

0.107 



0.010 



0.080 



0.016 



s(v E ) 

0.006 



0.010 



0.011 



0.005 




X]dimetilcarbonato+ 

x,dimetilcarbonato+ 








x 2 octano 
(298.15) 19 ' 20 


x 2 decano 

(298.15) 1 ’' 20 








7=0 

1.860 

-0.393 

0.981 

2.414 

-1.659 

1.749 








7.166x 

4.869x 

-1.486 

-7.869 

-9.628 

1.239x 







7=1 

10 4 

10' 4 

xlO' 3 

xlO' 6 

xlO 4 

10' 3 








1.008x 

1.933 x 

-1.506 

-4.756 

-1.137 

-2.153 







7=2 

10' 5 

10' 6 

xlO' 5 

xlO' 7 

xl0‘ 5 

xlO’ 6 








-8.121 

-4.869 

1.419x 

1.053 x 

1.849x 

4.73 lx 







7=3 

xl0‘ 9 

xlO' 8 

10' 7 

10' 7 

ÍO’ 8 

10' 8 







m 

0.001 


k = 
0.956 

0.002 


k = 
1.081 







stXi) 

0.005 



0.010 









s(v E ) 

0.005 



0.010 
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Tabla 2 (c). Valores de k y B i} obtenidos para cl modelo (29), en la correlación conjunta de los LLV isotérmico \ los v* a la 
temperatura indicada (T), de diferentes mezclas binarias, siendo s(\) la de$\iación estándar de la cantidad v. Aquí 0=Hg h RT) > k t se 
obtiene por optimización de (23) en el procedimiento de ajuste y k^(p) se calcula de (A.1.4), 




Bu 


floi - 


_ 


Bu 

B, ¡ 

B^ 

*■■ 

B : 


xietanol+.v : he\ano 

x,etanol+.r : hexano 

,rietanol+.r : heptano 

.r^etanol+r-heptano 

(313.15)'” 



(298.15) M 


(318.15) 


(298.15) 6 8 

7=0 

4.343 

-6.475 

3.655 

4.980 

-7.780 

4.076 

4.512 

-6.416 

3.951 

4.443 

-6.647 

3.839 


4.510x 

ío- 4 

5.109x 

8.450x 

3.576x 

-1.738 

8.450x 

-4.768 

-7.707 

-1.649 

4.711 x 

-5.746 

1.525 < 

7=1 

10* 3 

I0* 3 

10* 3 

xlO* 3 

10* 5 

xlO* 5 

xlO* 3 

xlO* : 

10* 5 

xlO 4 

10* 5 


-7.219 

1.633 x 

1.621 x 

-1.688 

-2.244 

1.562x 

-1.673 

-1.519 

-4.576 

-8.947 

5.725x 

1.527* 

7=2 

*10' 5 

10* 4 

10* 4 

xlO* 6 

xlO* 5 

10* 4 

xlO* 4 

xlO 4 

xlO* 5 

xlO* 5 

10* 5 

10 4 


5.178 x 

-1.365 

-1.253 

1.070x 

2.530x 

-1.177 

1.322x 

1.104 x 

9.452 - 

6.604 x 

2.393 x 

-9.442 

7=3 

10* 7 

xlO* 6 

xlO* 6 

10* 6 

10® 

xlO* 6 

10* 6 

ÍO 4 ’ 

10* 7 

10* 7 

10* 7 

xlO* 7 

s(0) 

0.021 


0>o>=1.767 

0.035 


k t 

(po)=\J61 

0.034 


k g 

(po)= 1.973 

0.014 


k t 

(p«)= 1.973 


0.707 


A\( 

2.074 



0.567 


M 

0.4S1 



s(yd 



/>o>=2.242 



/>0 >=2.256 



/>o >=2.514 



/>o >=2.519 

5(V F ) 

0.036 



0.057 



0.041 



0.042 




.rjpropan-1 -ol+x : heptano 

r, propan-1 -ol+.r : nonano 

r, butan-1 -ol+.r : hexano 

r,butan-l-ohr 

>heptano 



(298.15) 11 '* 


(298.15) ■ J 


(298.15) ' 5 


(313.15P 16 

7=0 

3.827 

-5.018 

2.911 

4.048 

-6.168 

3.739 

3.075 

-4.049 

2.608 

3.687 

-5.486 

2.8S1 


-6.413 

-8.451 

-1.211 

8.320x 

10* 5 

8.452x 

1.545 x 

-1.806 

-1.892 

1.152x 

2.111 x 

7.761 x 

8.314* 

7=1 

xlO* 5 

xlO* 5 

xlO* 3 

10* 3 

10* : 

xlO* 3 

xlO* 3 

10* 3 

10* 5 

10* 3 

I0* : 


-1.508 

-1.562 

-5.966 

9.552x 

1.050x 

1.513x 

-1.630 

-1.6S8 

1.619x 

-1.257 

-1.49 7 

-1.481 

7=2 

xlO* 4 

xlO* 4 

xlO* 5 

10* 5 

10* 4 

10* 4 

xlO* 4 

xlO 4 

10 4 

xlO 4 

xlO 4 

xlO 4 


1.245x 

1.207x 

4.953 x 

-8.420 

-1.097 

-1.398 

1.156x 

l.lOOx 

-1.293 

S.337x 

8.096x 

-1.652 

7=3 

10 -6 

10* 6 

10* 7 

x|0* 7 

xlO* 6 

xlO- 6 

10* 6 

10* 6 

xlO* 6 

10* 7 

10' 7 

xlO* 6 

s(Q) 

0.014 


k t 

(/zoH 1 -615 

0.016 


(/?o)= 1.967 

0.008 


(po)=1.219 

0.030 


k t 

(poHl .370 


0.266 


u 

0.199 


k,( 

0.322 


<\( 

0.754 


u 

s(y t ) 



po>=1.963 



/>o >=2.391 



/>o>= 1431 



/>oFl.612 

S^) 

0.021 



0.033 



0.016 



0.026 




Xiisobutil alcohol+.x : he.\ano 

.r|isobutil alcohol+r : heptano 

r,isobutil alcohol+r : heptano 

r, d imeti lcarbonato-r.whexano 



(298.15) 1417 


(298.15) 

2.1S 


(313.15) 6 

(298.15) l, -’°" 

7=0 

1.573 

-0.571 

0.162 

3.272 

-3.711 

2.142 

3.488 

-4.666 

2.234 

2.526 

-1.736 

0.931 


7.S46x 

-4.960 

-1.058 

-7.266 

-8.323 

-4.403 

-2.369 

2.523x 

6.130x 

3.422x 

-4.601 

6.325 x 


xlO* 5 

xlO* : 

xlO* 5 

xlO* 5 

xlO* 5 

xlO* 5 

10* 5 

10* : 

10* 5 

xlO* 5 

10* 5 

7=1 

10* :> 

6.624 x 
10* 5 

-1.470 

-1.409 

-1.559 

-1.567 

-9.109 

-1.686 

4.938x 

-7.921 

2.197x 

-7.221 

S. 152 x 

7=2 

xlO* 5 

xlO* 4 

xlO* 4 

xlO* 4 

xlO* 5 

xlO* 4 

10* 5 

xlO* 5 

10* 5 

xlO* 5 

10* 5 


-6.559 

1.778x 

10* 7 

1.302x 

1.309x 

10* 6 

1.205x 

8.062 x 

1.243x 

-5.177 

-1.402 

-3.890 

4.255 

-6.228 

7=3 

xlO* 7 

10* 6 

10* 6 

10* 7 

10 -6 

xlO* 7 

xlO* 6 

xlO* 8 

xlO* 7 

xlO*' 

^(0 

0.007 








k t 



k t 



(P»)= 1.212 

0.008 


(/?o)= 1.358 

0.025 


(/>o >=1.361 

0.002 


(po)= 1.264 


0.097 



0.103 


*>< 

0.384 


M 

0.027 


M 

5(^0 



PoH-416 



/?o>= 1.588 



/>o>=l -594 



/>o>=l .554 

5 (v r ) 

0.010 



0.020 



0.022 



0.021 




x,dimetilcarbonato+.r-octano 

.r^dimetilcarbonato+r-decano 








(298.I5) 193 ’ 

(298.15) l, ’°’ 







7=0 

3.018 

-3.137 

1.569 

3.817 

-5.738 

3.543 








-4.675 

-1.989 

1.585 x 

4.225 x 

8.265x 

2.277x 








xlO* 4 

xlO* 5 

10* : 

10* : 







7=1 


10’ 

10' 








-2.195 

-6.273 

1.248x 

1.3 55 x 

1.164 x 
10 4 

1.265 x 







7=2 

xlO* 5 

xlO* 5 

10* 4 

10 4 

10* 4 








3.095 x 

2.181 

-1.189 

-8.555 

-1.356 

-1.361 







7=3 

10* 7 

xlO* 7 

xlO* 6 

xlO* 7 

xlO* 6 

xlO* 6 







*(0 

0.006 


(p»)= 1.720 

0.027 


(A.H.163 








0.057 



0.165 


<\( 







5(/¡) 



/?oH.931 



/>o>=2.313 







5(V h ) 

0.024 



0.045 
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4. Conclusiones 

En este trabajo se presenta un modelo matemático flexible para la función de Gibbs 
adimensional g h /RT=y/ { ( x t p ) del que surgen otras expresiones derivadas para el volumen de exceso 
v E y para la función de Helmholtz a E de soluciones. Se pretende complementar un planteamiento 
realizado en trabajos anteriores [1-3] donde se desarrolló y analizó la aplicación de una expresión para la 

función de Gibbs con la temperatura, g E /RT—y/ 2 {x t T ). No debe olvidarse que ambas funciones se 
corresponderían con las parciales del objetivo último, la consecución de la relación 
g E /RT=y/(x t p,T ). Las relaciones que fundamentan este trabajo, indicadas en forma genérica con las 
ecuaciones (19) a (22), permiten establecer varias posibilidades de utilización en base al número de 
parámetros ajustables y también, al considerar o no, la dependencia con la presión del parámetro h X] , 
parámetro que caracteriza a la fracción activa “z” a través de la propiedad que se estudia. 



Figura 2 (a,b). Representación de valores experimentales, refs. [4,7] y curvas de ajuste obtenidas 
con el modelo simplifícado (26) para la mezcla (rietanol+j^hexano) a 318.15 K, con los parámetros 
presentados en la Tabla l(b). (a): (O), g E /RT vs X\; (A), y vs xú (b): valores de g E (V) y v E (O) y 
curvas de ajuste (-); (-), valores de a E , ecuación (28), y (-) h E de ref. [4]. 


Para validar el modelo propuesto se han escogido un número de sistemas binarios para los que la 
literatura presenta datos de ELV isotérmicos y volúmenes de exceso, en algunos casos en condiciones 
diferentes entre ellos, de presión y temperatura. La utilización del modelo aquí propuesto da lugar a seis 
casos posibles cuyos resultados se recogen en las tablas l(a-c) y 2(a-c). E1 resumen de las diferentes 
posibilidades estable ce que el empleo de un modelo extendido debería ser útil solo para ciertos casos 
califkados de “difíciles”, ya que otra cuestión es que la dependencia de h Xi con la presión no sea 
significativa (ya que su consideración ha sido en base a un acercamiento real del problema). Es claro que 
si el cálculo del parámetro k\¡ se realiza por optimización como un coeficiente más del modelo en el 
procedimiento global de ajuste, surge una mejora de los resultados de la correlación. Por ello, al menos 
por los resultados que aquí se obtienen, la propuesta concreta más simple sería utilizar la ecuación (30) 
como principal, haciendo que k v =h B sean sometidos al proceso de correlación con el AG. E1 caso de 
reducción de parámetros quedará a criterio del investigador y su empleo parece obvio. 

En la figura 2(a,b) se muestra gráficamente uno de los ejemplos utilizados, el del sistema binario 
(etanol+hexano) a 318,15 K, con el modelo simplificado, remarcado en las tablas l(a-c) y 2(a-c) 
observándose en las fíguras la bondad cualitativa del ajuste de datos experimentales y la validez del 
modelo y procedimiento empleado, tanto para las curvas de equilibrio líquido-vapor como para la de v E . 
En la fígura 2(b) se aprecia una superposición de los datos de g E y a ] , ya que los valores del término (pv ) 
resultan despreciables frente a g E . 
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6. Apéndice 

6. 1. Cálculo de los parámetros k X] 

Para la correlación de volúmenes de exceso se propuso el cálculo de los coeficientes k u de (9) a 
través del cociente de volúmenes molares de las sustancias puras y en las mismas condiciones de trabajo, 
esto es: 


v;(t.p) 

(*) 

’pJT.pí 

V°(T.p) 

Wj 

yjT. P )_ 


este cociente se identificó en anteriores artículos [1-4] como y para ellos se estableció la constancia con 
la presión. La variación de dicho parámetro con la temperatura se planteó en [1] utilizando el coeficiente 
de expansión térmica a de las sustancias puras. La expresión final para A, se consiguió introduciendo los 
valores de estos coeficientes en la (A.l) y fúe la siguiente: 


k,(T) = k,(T) = ^ = 


M i P\ 


K M\ P° 


M, 


M, 


P\ -e 


-a t (T-T 0 ) f 


¡ y 


o -ayT-T 0 ) 


P, • 


V 


M. 


P\ 


Kp 


(a-a,\T-T 0 ) _ 


j / 


(A.2) 


= k v (T 0 )e <artt ' KT ~ T ° ) 

siendo M x , M } , p° x , , respecti\ amente, las masas moleculares para las sustancias i y j, y las 

correspondientes densidades a una temperatura de referencia, T 0 . Realizando un estudio paralelo, la 
dependencia de los volúmenes con la presión puede llevarse a cabo mediante la introducción del 
coefíciente de compresibilidad isotermo, parámetro característico de las sustancias y establecido 
conceptualmente de la siguiente manera. 


1 

'dv) 


f dln p') 

V 


T 

l cp J 
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y cuya integración da lugar a: p = p 0 e^ (p Po) . La introducción de esta relación en la (A.l) para 
dos sustancias i y j daría lugar a la siguiente expresión para k v como función de la presión. E1 valor de k h 
se obtiene a partir de la relación establecida en un trabajo anterior [3]. 



( M i) 

p° . e ^p-p»> 

í 


V° M, p. 

UJ 

° Ps(p-Po) 
p. e 


U°J 


¿P-P-^p-Po) 


(A.4) 


= K(P a )e P¡ 


(ft-PiXp-Po) 


siendo ahora p° y p° } , las correspondientes densidades a una presión de referencia, p 0 . Ambos 

coeficientes, el de expansión ténnica isobárico a , y el de compresibilidad isotermo fi\ dependen de la 
temperatura y de la presión, por lo que, las expresiones correspondientes (A.2) y (A.4) podrían ser aún 
más complejas; sin embargo, la consideración de un valor promedio para ambos coeficientes, ay /?, en un 
amplio interv alo de temperaturas y presiones puede ser aceptado. 

Analizada la dependencia de la fracción activa “z” con la temperatura [3], comprobaremos la 
influencia de “z” con la presión, introduciendo la derivada matemática del parámetro k x con relación a p 
que puede obtenerse de manera sencilla a partir de (A.4) las cuales pueden sustituir a la correspondiente 
cantidad en las ecuaciones (19), (21) y (22) para completar el modelo. 

Como aspecto práctico comentar que los coeficientes de compresibilidad isotermo son valores de 
muy pequeño orden, por lo que el exponente de (A.4) tiende a ser cero y la expresión final para k v es 
equivalente a la que se obtendría sin considerar la variación con p. Sin embargo, la idea de este trabajo es 
plantear un estudio riguroso del problema con el fin de generalizar su aplicación teniendo en cuenta los 
aspectos formales del estudio termodinámico del equilibrio entre fases. 
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Figura A.l. Variación de los volúmenes 
molares de sustancias puras, alcanoles y 
alcanos, con la presión 


Figura A.2. Representación del parámetro k x (p) 
vs p para varias mezclas: C u H 2 u+i(OH)+C n H 2 n+ 2 , 
indicándose como u+n. 
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La aplicación numérica de este apartado se ha realizado con datos de densidad de los primeros 
alcanoles e hidrocarburos, en función de la densidad, obtenidos de la base de datos del TRC [21]. Así, en 
la figura (A.l) se recoge la variación de los volúmenes molares con la presión de un conjunto de 
sustancias puras pertenecientes a las mismas familias de los compuestos utilizados como ejemplo en este 
trabajo. En esta gráfica se aprecia claramente una ligera pendiente negativa, a altos valores de presión, de 
la distribución de los K°(p). Sin embargo, dicha pendiente se hace casi nula cuando se evalúan los 
cocientes, representados en la figura (A.2), manteniéndose prácticamente constante los valores de k v con 
la presión para una determinada mezcla binaria. A pesar de ello, en esta serie de trabajos en la que se 
intenta plantear un procedimiento riguroso para el tratamiento de las propiedades termodinámicas de 
sistemas en equilibrio, se ha considerado la correspondiente variación de dicho parámetro con la relación 
(A.4). 

Evidentemente, los cambios correspondientes del parámetro entálpico k h con la presión serían 
introducidos de igual forma a como se detalló en un trabajo anterior [3], a través de Recordemos 
nuevamente esta dependencia, 


K = 


í Y V 3 

?i 


k; 3 (T) 


expresión que también podría haberse escrito como 


= 


V- \ 2n 

<h_ 'j 


V r ¡ J 


k; 3 ( P ,T)=k h {Po,T 0 )k; l3 {p,T) 


(A.5) 


(A.6) 


siendo la relación establecida en (A.2) y/o en (A.4) la que podría utilizarse para plantear, respecti\ amente, 
la variación de k x con la temperatura o con la presión. Basándonos en la relación (A.6) y las 
correspondientes para k vy proponemos aquí que la \ariación del parámetro A g se represente por una 
ecuación análoga, es decir, 

k % =k g (Po j. Xv 3 (P ’ T ) (A.7) 


lo que sugiere que podría iniciarse el cálculo con un \ alor en determinadas condiciones de presión y 
temperatura y hacer variar dicho parámetro, bien con la presión o con la temperatura, de acuerdo a la 
(A.7). 
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